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Abstract 

We show that the Collins-Soper-Sterman resummation approach to the 
derivation of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi equation is gauge 
invariant. The special gauge- dependent parton distribution function em- 
ployed in the resummation technique is expressed as the convolution of an 
infrared finite function with the standard distribution function. By means 
of this convolution relation, we explain how the technique works in summing 
large logarithmic corrections, and how the gauge invariance of the special 
distribution function is restored. 
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1. Introduction 



It is known that a quark distribution function for a hadron in the minimal 
subtraction scheme is defined by 

0(£,/Vm) = / ^e^ + y-{p\q(y-)^ + Pe i Sf dzv ' A ^\(0)\^ , (1) 

where 7 + is a Dirac matrix, and \p) denotes the incoming hadron with 
the momentum p M = p + 5 t±+ . Averages over spins and colors are under- 
stood. The above expression, with the presence of the path-ordered expo- 
nential Pe l f dzv ' A ^ zv \ = being a light-like vector, is gauge invariant. 
Through this exponential, <fi collects the collinear divergences of radiative 
correcctions to a QCD process. While the infrared finite piece of radiative 
corrections is absorbed into a hard scattering amplitude. The factorization 
formula for a cross section (or a structure function) is then expressed as the 
convolution of these two factors. ///ra) describes the probability that a 
parton carries the fractional momentum £p at the factorization (or renormal- 
ization) scale /j,, with m an infrared cutoff for the collinear divergences, such 
as quark mass. The argument fi/m denotes the large logarithms ln(/i/m) 
from the collinear divergences, which will be summed by the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation 

To organize the large logarithms contained in the distribution function 
using the Collins-Soper-Sterman resummation technique 0, the vector v' is 
replaced by an arbitrary vector n, n 2 ^ 0, leading to a n-dependent distri- 
bution function, 

z///i, /*/"») = j^e- iip+v ~{p\q{y-)\l + Pe i! ° ^ (jB) g(0)|p> ■ (2) 

We have made explicit the additional argument v = J {p ■ n) 2 /\n 2 \ of (p^ n \ 
which appears as a ratio because of the scale invariance of Eq. (0) in n 
Unfortunately, this replacement spoils the gauge invariance of Eq. ([!]), since, 
with n being arbitrary, the end point of the path, on which the gauge field 
A is evaluated, does not coincide with the coordinate y~ of the quark field q. 
It will be shown that the scale v serves as an ultraviolet cutoff for the loop 
integrals associated with the collinear gluons, so that (j) and (j)^ possess the 
same infrared structures but different ultraviolet structures. 
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We shall demonstrate that the evolution kernel of the DGLAP equation, 
obtained from the resummation technique, turns out to be n- independent, 
i.e., gauge invariant. After deriving the equation, n is brought back to the 
light cone. In the n — > v' limit v diverges, and ln(z///i) corresponds to an 
ultraviolet pole in dimensional regularization. This pole is removed in a 
renormalization scheme, and <^ n > reduces to the gauge invariant distribution 
function <p. We also show that (p^ approaches <p as v = fi, for which the 
extra logarithms ln(z// fx) vanish. That is, the arbitrary vector n appears only 
at the intermediate stage of the formalism, and as an auxiliary tool. The 
gauge invariance of the resummation approach to the evolution equation is 
thus guaranteed. To make our reasoning concrete, we derive the DGLAP 
equation in both covariant and axial gauges. The above discussion can be 
extended to the gluon distribution function directly, whose evolution in the 
small momentum fraction is governed by the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation [[§. 

We explore the relation between the working definition (p^ and the stan- 
dard definition <ft: the former is expressed as the convolution of an infrared 
finite //-dependent function with the latter, because they contain the same 
nonperturbative information. Hence, varying n can be understood as dif- 
ferent ways to partition the infrared finite contribution to a cross section 
between the parton distribution function and the hard scattering amplitude. 
The n dependence is then in fact a factorization scheme dependence. By 
means of this convolution relation, we analyze in details how the resumma- 
tion for 0( n ) is accomplished, and how becomes <p in the limit n — > v'. 

2. The DGLAP Equation 

Take deep inelastic scattering of a hadron as an example. We first explain 
how to factorize collinear gluons into the quark distribution function, and 
justify the replacement of v' by n in the collinear kinematic region ||. It 
is known that the collinear region with loop momenta of radiative gluons 
parallel to p is important, from which large logarithms arise 0. In this 
region the partial integrand for a loop diagram, with the radiative gluon 
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attaching the scattered quark line, is approximated by 



where the scattered quark momentum p' possesses a large minus component 
p ~ in the large x limit. The preceding factor jp' comes from the final-state 
cut. Equation fl3|) indicates that the scattered quark line, the collinear gluons 
attach, can be replaced by an eikonal line in the direction v'\ 1/(V • I) is the 
Feynman rule for an eikonal propagator, and v' for a vertex on the eikonal 
line, / being the momentum flowing through it. It is straightforward to show 
that these rules are exactly produced by the path-ordered exponential in 
Eq. ([!]). With the eikonalization, the collinear gluons are factorized into the 
quark distribution funtion cf). 

We may employ the further approximation, 



v' ■ I l + n /+ n /+ + n + l — • n • I 

The smaller components l~ and It have been added to the denominator, and 
the terms proportional to the components n + and tit of an arbitrary vector 
n, which give vanishing contributions when contracted with a vertex in the 
distribution function, have been included in the numerator. Similarly, the 
path-ordered exponential in generates the Feynman rules n^/ (n ■ I) in 
Eq. (U). Because (j)^ and <p coincide with each other in the collinear region, 
they contain the same nonperturbative information, and their difference is 
infrared finite as stated in the Introduction. A general diagram for (p^ is 
exhibited in Fig. 1(a), where the exponential is represented by an eikonal line 
along n with collinear gluons attaching it. 

In the resummation framework the DGLAP equation can be obtained by 
studying the derivative p + dcfi^ / dp + in the covariant gauge d ■ A = 0. Due 
to the argument u, we have the relation 0. [| . 

P + ^ n) = -—v a -^-cpM • (5) 
dp + v ■ n dn a 

with = a vector along p. Since p flows through all the quark and gluon 
lines in the distribution function, while n appears only in the exponential, 
Eq. (Kf) simplifies the analysis of the p + dependence of (f)( n \ This is the 
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motivation to introduce the arbitrary vector n. The differentiation of Eq. 
with respect to n a , 



n 



d n. 



v ■ n dn n n ■ I 



n 2 (vl 
v ■ n \n ■ I 

n ■ I ' 



n ■ I 



gives 



d 
dp^ 



(n) 



Un) 



(6) 



(7) 



which is described by Fig. 2(a). A summation over different attachments of 
the symbol x, which represents the special vertex n M defined by the second 
expression in Eq. ([]), is understood. The coefficient 2 comes from the equality 
of the new functions (f)^ with the special vertex on either side of the final- 
state cut. 

If the loop momentum I is parallel to p, the factor v ■ I vanishes, and n M is 
proportional to v^. When this is contracted with a vertex in (f)( n \ where 
all momenta are mainly parallel to p, the resultant contribution diminishes. 
Therefore, the leading regions of I are soft and hard, in which the subdiagram 
containing the special vertex can be factorized, leading to 



(8) 

with the functions K and G absorbing the soft and ultraviolet divergences, 
respectively. The 0(a s ) contributions to K from Fig. 2(b), where the eikonal 
approximation for the valence quark propagator has been made, and to G 
from Fig. 2(c), where the soft subtraction ensures a hard loop momentum 
flow, are written as 



K 

G 
C F 



d 4 ~H 



In e f « 

19 C F n J — 



(27r) 4 ~ e n-lv-l 

d A ~n 



X 



n,. 



(27r) 4 " e n ■ IP 
4/3 being a color factor. 



/ 



(p + : 



+ 2Tti8(l 2 )5 \Z-x 



p+ 



vl 



(9) 
(10) 
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A straightforward calculation gives 

7/ 

(11) 



' + In -5(1 - x/i) 



(1-^/0+ H 



K = ^>C, 
G = -^lc F ln^5(l-x/0 , (12) 

where constants of order unity have been dropped. In the considered region 
with x — > 1 the logarithm ln(£z///i) in G can be replaced by ln(i///x). Adding 
the above expressions, we have 

K(x/t, u/fi) + G{x/t, u/im) = ^rC F 1 , (13) 

^ (1 - + 

where the argument of a s has been set to the characteristic scale v of the 
evolution kernel K + G. It is observed that the gauge factors u have cancelled 
between K and G. 

Inserting Eq. ( |13"D into @, Eq. (^) becomes 



dp+ ' 7T A £ (1 ~X/Cl + 

(14) 

The solution to the above differential equation can be written as 

{n) (x,l,/i/m) = <f^ n \x,A/fi,fi/m) 

+ r?<2) Cr /** *^w ({ifl//Ii/l/m) , 

JA fl TV Jx i (1-X/0 + 

(15) 

with A an arbitrary scale. We shall explain in Sect. 3. that (p^ coincides 
with <p as v — i.e., as the logarithms \n(i//fi) vanish: 

(n) (x, 1, n/m) = (f>(x, n/m) . (16) 

Differentiating Eq. ([15]) with respect to /i, and employing the RG equation 
shown in Sect. 3, 

pfy ln) M£)> A Oi)/t*,t*/m) = -2A 9 ^)(x(0,A(/i)/^/i/m) , (17) 



with X q = —a s jn the quark anomalous dimension, we obtain 



jjL—<p{x,jj,/m) = Lp 

d\i n 



-<!>(£, /V m ) ~ 2A g (/x)0(x, n/ m) 



(18) 

At this step, the n dependence disappears completely. The above formula is 
then identified as the DGLAP equation, 



fi-^-(j)(x, fi/m) 
dfj, 



0£ a {p) f 1 d£ 



7T 



(19) 



where the evolution kernel P is given by 



P(x)=C, 



1-x) 



+ ^(l~x) 



(20) 



Next we derive the DGLAP equation in the axial gauge || ||, in which 
the relevant definitions for the quark distribution function are 



0M(£,/i/m) 
W(£,iV/i,/i/m) 



^-e-^ + y-(p\q(y-)^ + q(0)\p)\ v , A=0 , (21) 
^ e -^ + y-{ P \q(y-)^+q(0)\p)\n-A=o , (22) 



as described by Fig. 1(b). The path-ordered exponential is equal to the 
identity in this gauge, and eikonal lines collecting the collinear gluons are 
absent. Similarly, we work on (j)^ n > when performing the resummation. The 
n dependence goes into the gluon propagator, (— i/l 2 )N ilv {l) y with 



n ■ I [n ■ iy 



(23) 



Because of the scale invariance of in n as indicated by Eq. (p3|) , <j>^ 
depends on the scale u, and thus Eq. @ holds. The operator d/dn a , now 
applying to the gluon propagator, gives 



71 -v*N>" 



n 2 v, 



v ■ n dn r 



— (FN au + l u N» a ) 



v ■ nn ■ I 

v a {FN av + l v N m ) . 



(24) 
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The loop momentum (l u ) carried by the differentiated gluon contracts with 
the vertex the gluon attaches, which is then replaced by the special vertex 
v a defined by the second line of Eq. ( P^ ) . 

The contraction of l v hints the application of the Ward identities || §], 
such as 



(k + iy 

for the quark-gluon vertex, and 



i(/k+ /) f i _i ]k i(/k+ /) 



2 ' 



-iN a "(k + l) 

(fc + 2 



-iN x ^(k) 



k 2 



(k + l) 



-iN a "f(k) -iN a "f(k + l) 



k 2 



(k + lf 



(25) 



(26) 



for the triple-gluon vertex V ^ v \. Similar identities for other types of vertices 
can be derived easily. Summing all the diagrams with different differentiated 
gluons, those embedding the special vertices cancel by pairs, leaving the 
one with the special vertex moving to the outer end of the quark line 0. 
We obtain the same formula in Eq. (0), described by Fig. 3(a), where the 
new function ^ n > contains one special vertex represented by a square. The 
coefficient 2 comes from the equality of the new functions with the special 
vertex on either of the two quark lines. 

The important regions of the loop momentum / flowing through the spe- 
cial vertex are also soft and hard, since the vector n does not lie on the 
light cone, and the collinear enhancements are suppressed. Similarly, in the 
above leading regions can be factorized into the convolution of the sub- 
diagram containing the special vertex with the original distribution function 
(fft 1 ' , leading to Eq. (^). In this case the functions K and G, corresponding 
to Figs. 3(b) and 3(c), are written as 



K 



■ In e 

ig Cf(J, 



d 4 ~n 



(27r) 4 " e v ■ I 



- x) 
i 2 



G = -ig 2 C F n e J 



(2tt 



|4-£ 



v-l 



I 2 



(27) 
(28) 



It is easy to show that Eqs. (|27j ) and ( p8|) reduce to Eqs. ( ITT] ) and (|T2|), 
respectively. Then following the steps from Eq. (|13"D to Eq. (0), we derive 
the same DGLAP equation in the axial gauge. This is natural, because 
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the evolution of a parton distribution function is measurable, and should be 
independent of the gauge we adopt. 



3. Relation of dM> and 



In this section we explore the relation between the working definition (j)^ n > 
of the quark distribution function employed in the resummation technique 
and the standard definition <p. As explained in Sect. 2, (p^ and <\> contain 
the same collinear divergences from gluon momenta parallel to the hadron 
momentem, and thus the same nonperturbative information. Take the 0(a s ) 
correction with one end of a real gluon attaching the quark line and the other 
end attaching the eikonal line as an example. The loop integral associated 
with <f) is written as 

t 2 r e f d 4 ~H (x fi+ Jt) j6' n fl ( l+\ 

1 = 9 J (27T)4-e (xp + I) V ■ l 27T6{l } H ' 



a s /i x 



d^l Tir -— 2 . (29) 



7i (2vr) 1 - e 2(l -x) J U^ + m 2 

where m is the infrared cutoff mentioned before, and /x is the renormalization 
scale. The corresponding integral I^ n > associated with (j)^ n > is obtained by 
replacing v' in Eq. ( P^ ) by n, which is given by, 

An) [jtj 2x{\-x)v 2 

~ 2tx 2 J (/ 2 + m2) [Z 2 + 4(1 _ x y v 2} ■ ^) 

The collinear divergences come from the region with small It, in which 1^ 
reduces to I as indicated by Eqs. fl2~9|) and (p0|). However, / is ultraviolet 
divergent, while is not, because of the extra denominator l T + 4(1 — 
x) 2 v 2 . Therefore, (j)^ and </> possess the same collinear structure but different 
ultraviolet structure. 

Performing the integrations over It, we obtain 

/ = ^-(1+lnii) , (31) 
2tt 1 -x V2e mj y ' 

jin) = 21-E-] n IL } (32) 

2tt 1 — x m 
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where constants of order unity have been dropped. It is found that I contains 
a logarithm ln(/x/m), and I^ 1 ' contains ln(z//m) with the same coefficient. 
Corrections in 0( n ) and without gluons attaching the eikonal line, such as 
the self-energy corrections, give the same logarithms ln(/x/m). By splitting 
the logarithm in I^ 1 into 

ln(z//m) = ln(z///i) + ln(/i/m) , (33) 

it is easy to understand why (j)^ depends on the argument u/fM due to the 
replacement of v' by n. The resummation technique, summing the first loga- 
rithm in the above expression to all orders, leads to Eq. (|T4]) in the same form 
as the (partial) DGLAP equation, that sums the second logarithm, since the 
two terms possess the same coefficients. The choice v = /j,, diminishing the 
first term, renders ( - T1 ) coincide with as in Eq. flTBl) , because the second 
term, combined with those from other corrections without gluons attaching 
the eikonal line, gives the complete logarithms ln(///m) in 0. It is also easy to 
realize why the RG method, applying to 0(™)(x(£), A (//)///, n/m), sums only 
ln(/i/m) from other corrections as in Eq. (|i~7|), since the corrections with glu- 
ons attaching the eikonal line, proportional to ln(z//m), are /i- independent. 
At last, in the limit n — *■ v' the cutoff v diverges, and ln(z///i) corresponds to 
the ultraviolet pole in /. This pole is removed by renormalization, and (p^ 
approaches 0: 

lim ( - n * ) (x,u/ fi, fi/m) = 0(x, /i/m) . (34) 

n —*v' 

The infrared finite difference between <n - ) and 0, indicated by 1^ — I oc 
ln(z//yu) from Eqs. (|3T|) and fl32|), can be computed perturbatively. In terms 



of Feynman diagrams in the covariant gauge, it is the difference between 
the diagrams with radiative gluons attaching the eikonal line along n and 
the diagrams with gluons attaching the eikonal line along v'. The above 
argument leads to the convolution formula, 

0^(x,z// yU , yU /m) = £ ^DW^v/M&p/m) , (35) 

where the infrared finite piece D, independent of the infrared regulator m, 
denotes the difference mentioned above. The graphic defintion of D is given 
in Fig. 4, where the denominator, with the eikonal line along v', subtracts 
the collinear divergences involved in the numerator with the eikonal line 
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along n. All the diagrams without gluons attaching the eikonal lines cancel 
between the numerator and the denominator. The bubbles contain infinite 
many gluon exchanges, and the external lines represent quarks. Because of 
D, 0( n ) absorbes additional finite contributions compared to <p. Hence, they 
can be regarded as the distribution functions defined in different factorization 
schemes. Varying the vector n then means varying the factorization scheme. 
Based on Eq. fl35f ), the differentiation d/dp + or d/dn in fact applies to D: 



p + -^<j> in \x 7 is/v,fi/m) = I'lp+^^^^UH . (36) 
The similar resummation procedures give 

(37) 

Using the variable change z = y/£, and interchanging the integrations over 
£ and over y, the above equation is reexpressed as 

p +_^o.> (li „ /m) = £ *L J> |_l_ z , ( , /c , w , m< , /m) . 

(38) 

Employing the definition in Eq. (|36|), we arrive at Eq. fllij) . The remaining 
steps of the resummation and Eq. (31) in the n — > v' limit then lead to 



the DGLAP equation. The above formulas elucidate the statement that the 
vector n is an auxiliary tool, and appears only at the intermediate stage of 
the derivation. Hence, the gauge invariance of the resummation technique is 
guaranteed. 

The discussion presented here can be generalized to the resummation- 
based derivation of the BFKL equation for the gluon distribution function 
F(x, kr) easily, which describes the probability that a gluon from the hadron 
carries the longitudinal momemtum xp + and the transverse momentum kx- 
However, the relation between the working definition and the standard 
definition F is a /c^-factorization formula , instead of a collinear factoriza- 
tion formula in Eq. (p5|). The details will be published elsewhere. 



4. Conclusion 
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In this paper we have demonstrated how to derive the DGLAP equation 
using the Collins-Soper-Sterman resummation technique, and that the results 
are independent of the gauge we employed. Though we started with the n- 
dependent parton distribution function, the evolution kernel turns out to be 
gauge invariant. We have also explored the relation between the n-dependent 
and standard distribution functions. From this relation, we explained why 
the resummation technique is a successful method to the summation of large 
logarithms, and showed that the n-dependent definition reduces to the stan- 
dard one by either taking the n — ► v' limit or setting the gauge factor v to the 
renormalization scale u. This work provides a solid ground for our previous 
studies on the resummation approach to evolution equations || |(| . 

This work was supported by the National Science Council of R.O.C. under 
the Grant No. NSC-87-2112-M-006-018. 
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Figure Captions 

FIG. 1. Definition of the quark distribution function in (a) the covariant 
gauge and in (b) the axial gauge. 

FIG. 2. (a) The derivative p + d(J)^/dp + in the covariant gauge, (b) The 
0(a s ) function K. (c) The 0(a s ) function G. 

FIG. 3. (a) The derivative p + d(j)^/dp + in the axial gauge, (b) The 0(a s ) 
function K. (c) The 0(a s ) function G. 

FIG. 4. The graphic definition of the function D. 
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